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Abstract 



We consider partitions of the positive integer n whose parts satisfy 
the following condition: for a given sequence of non-negative numbers 
2 ' {bk}k>i, a part of size k appears in exactly 6^ possible types. Assuming 

that a weighted partition is selected uniformly at random from the set of 
all such partitions, we study the asymptotic behavior of the largest part 
X„. Let D{s) = X^^j^ bkk'" , s — a + iy, be the Dirichlet generating series 
of the weights bk- Under certain fairly general assumptions Meinardus 
(1954) has obtained the asymptotic of the total number of such partitions 
as n — >■ cxD. Using Meinardus scheme of conditions, we prove that X„ 
i~^ , appropriately normalized, converges weakly to a random variable having 

^^ ' Gumbel's distribution (i.e. its distribution function equals e^*^ , — cxj < 

l»^ , t < oo). This limit theorem extends some known results on particular 

f- ^ ■ types of partitions and on the Bose-Einstein model of ideal gas. 

Mathematics Subject classifications: 05A17, 60C05, 60F05 
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H . 1 Introduction and Statement of the Result 

A weighted partition of the positive integer n is a multiset of size n whose de- 
composition into a union of disjoint components (parts) satisfies the following 
condition: for a given sequence of non-negative numbers {6fc}fc>i, a part of size 
k appears in exactly one of bk possible types. For more details on properties 
of multisets, we refer the reader e.g. to [3]. Weighted partitions are also asso- 
ciated with the generalized Bose-Einstein model of ideal gas, where n(= E) is 
interpreted as the total energy of the system of particles. The weights bk,k> 1, 
are viewed as counts of the distinct positions of the particles in the state space, 
where a particle in a given position has (rescaled) energy k (for more details on 
the relationship between combinatorial partitions and various models of ideal 



gas, see [19]). From combinatorial point of view, it is fairly natural to assume 
that bk,k > 1, are integers (see e.g. the "money changing problem" discussed in 
detail in [22; Sect. 3.15]). On the other hand, it turns out that this requirement 
is not necessary for the analytical approach used in this paper. That is why, we 
assume that bk,k > 1, are real non-negative numbers. 

For a given sequence b = {bk, k > 1}, let Vb{n) be the set of all weighted 
partitions of the positive integer n and let Pb{n) =\ Vb{n) \ be its cardinality. It 
is known that the generating function fb{x) of the numbers Pb(n) is of Euler's 
type, namely, 

oo oo 

fbix) = l + Y,Pb{n)x^ = l[(^-x'')-'\ \x\<l (1.1) 

»i=i fc=i 

(see [22; Sect. 3.14]). We introduce the uniform probability measure P — Pn.b 
on the set of weighted partitions of n assuming that the probability l/pb{n) 
is assigned to each n-partition with weight sequence b. In this paper we focus 
on the size of the largest part Xn of a random weighted partition of n. With 
respect to the probability measure P, Xn becomes a random variable, defined 
on the set Vb{n). It is also well-known that 

oo m 

/,„,6(a:) = l + ^P6(r^)P(X„<m)a;" = [|(l-a;'=)-^^ m>l (1.2) 

n=l k=l 

(see [22; Sect. 3.15]). 

The asymptotic behavior of the combinatorial numbers Pb{n) (the Taylor 
coefficients in (jl.ip ) will play important role in our further analysis. A fairly 
general scheme of assumptions on the parametric sequence b was proposed by 
Meinardus [11] (see also [2; Chap. 6]), who found an asymptotic expansion 
for the numbers Pbin) as n — >■ cx). His approach is based on considering two 
generating series: 

oo 

D{s) =Y,bkk-' , s = a + iy, (1.3) 

fc=i 
and 

oo 

G{z)=Y,bkz\ \z\<l. (1.4) 

fc=i 

Below we give Meinardus' scheme of conditions. Throughout the paper by '^{z) 
and '^{z) we denote the real and imaginary part of the complex number z, 
respectively. 

(Ml) The Dirichlet series (|1.3p converges in the half-plane a > p > Q and 
there is a constant Co G (0,1), such that the function D{s) has an analytic 
continuation to the half-plane {s : a > —Co} on which it is analytic except for 
the simple pole a.t s — p with residue A > 0. 

{M2) There exists a constant Ci > such that 

D{s)=0{\yf^), I y 1^00, 



uniformly for a > — Cq. 

(M3) There are constants e > and C2{— C-zi^) > 0), such that the function 
g{T) — G{e^'^), T = a + 2TTiu, u real and a > (see (|1.4p ) satisfies 

5R(5(t)) - 5(a) < -C2a-% | arg(T) |> ^4, ^| w |< 1/2, 

for enough small values of a. 

The first assumption (Ml) specifies the domain, say H, in which D{s) has an 
analytic continuation. The second is related to the asymptotic behavior of D(s), 
whenever |3(s)| — >■ 00. Functions, which are bounded by 0(|5(s)|'^), < c < 00, 
in certain domain, as |3(s)| ^ 00, are called functions of finite order. It is known 
that the sum of the Dirichlet series in (|1.3p satisfies the finite order property 
in its half-plane of convergence a > p (see e.g. [18; Sect. 9.4]). Meinardus 
second condition requires that the same holds for the analytic continuation of 
D(s) in the domain H. Finally, Meinardus third condition implies a bound on 
5ft(G(e^^)) (sec (|1.4p ) for certain specific complex values of r. In some cases its 
verification is technically complicated. Granovsky et al. [9] showed that it can 
be reformulated as follows: 

00 

V bke-''" sin^ {irku) > Csa"', < ^ <| m |< 1/2, 
fc=i 

for small enough a and some constants 6*2,6 > (C2 = ^2(6)). 

Moreover, they proved that this inequality holds for any sequence bk,k > 1, 
satisfying the inequality bk > Ck'^~^, k > ko, for some fco ^ 1 and CjV > 0. We 
notice that if 

bk - Ck'''\ fc > 1, (1.5) 

then D{s) = CC{s — i/ + 1), where C denotes the Riemann zeta function. There- 
fore, D{s) has a single pole at s = v with residue C > and a meromorphic 
analytic continuation to the whole complex plane [21; Sect. 13.13]. These facts 
show that Meinardus conditions (Mi) — (A/3) are satisfied by the weights (|1.5p 
with p — V and A = C . 

Throughout the paper we assume that conditions (Mi) — (M3)) are satisfied. 
Our aim is to determine asymptotically, as n — > cxd, the distribution of the 
maximal part size X„. Recalling ()1.2p . we also point out that our results may be 
interpreted in terms of the asymptotic of the combinatorial counts of partitions 
whose part sizes are < m, where the range of values of m is specified by the weak 
convergence of the random variable X„ to a non-degenerate random variable. In 
the brief review given below we summarize some known results on the limiting 
behavior of the random variable X„ . 

Consider first the classical case of linear integer partitions, where the weights 
satisfy 6fc = 1, fc > 1. This kind of partitions were broadly studied by many au- 
thors in many respects. Their graphical representations by Ferrers diagrams 
show that their total number of parts and their maximal part size Xn are iden- 
tically distributed for all n (see [2; Sect. 1.3]). Erdos and Lehner [6] were ap- 
parently the first who applied a probabilistic approach to the study of integer 



partitions. As a matter of fact, they found an appropriate normalization for X„ 
in this case and showed that 7rX„/(6n)^/^ — log ((6n)^/^/7r) converges weakly, as 
n — > cx), to a random variable having the extreme value (Gumbel's) distribution. 
The local version of their theorem was derived later by Auluck et al. [4] . Frist- 
edt [8] studied linear integer partitions using a transfer method to functionals of 
independent and geometrically distributed random variables. Among other re- 
sults, he obtained the limiting distribution of the fcth largest part size whenever 
k is fixed. Finally, we notice that among weighted integer partitions only the 
linear ones possess the property that number of parts and maximum part size 
are identically distributed. The limiting distribution of the number of parts in 
the general case of random weighted partitions under Meinardus scheme of con- 
ditions is studied in [13]. It turns out that the limiting distribution laws depend 
on particular ranges in which the parameter p varies (see condition (Mi)). 

Another important particular case of weighted partitions arises whenever 
bk — k,k > 1. It turns out that in this case the generating function fb{x) 
(see (II. ip ) enumerates the plane partitions. A plane partition of n > 1 is a 
matrix of non-negative integers arranged in non-increasing order from left to 
right and from top to bottom, so that their double sum of its elements equals 
n. Together with the largest part size X„, consider also the counts of the non- 
zero rows and columns of the matrix of a plane partition. It turns out that 
these three quantities measure the sizes of the corresponding solid diagram of a 
plane partition in the 3D space. (The solid diagram is a heap of n unit cubes 
placed in the first octant of a coordinate system in a SD space whose columns 
composed by stacked cubes have non-increasing heights along the x- and y- 
axis; the height of this heap along the z-axis is just Xn, the largest part size.) 
Similarly to Ferrers diagrams for linear integer partitions, the three sizes of this 
heap appear to be identically distributed for all n > 1 (for more details, see 
[17; p. 371]). Their joint limiting distribution was found in [15]. The marginal 
limiting distributions (including the limiting distribution of Xn) were obtained 
in [12]. For more details on various properties of plane partitions and their 
applications to combinatorics and analysis of algorithms, we refer the reader to 
[2; Chap. 11], [14; Chap. 11] and [17; Chap. 7]. 

Our study is also closely related to some recent results on the maximal 
particle energy in the Bose-Einstein model of ideal gas. The general setting 
and the probabilistic frame of problems from statistical mechanics and their 
relationship with enumerative combinatorics was given by Vershik [19]. In the 
context of the infinite product formula ()1.1|) , he studied the Bose-Einstein model 
by a family of probability measures fiv,v Cz (0, 1), defined on the set of all b- 
weighted partitions Vb = U„>o7^f)(n.). So, for a partition A = (Ai,...,Ai) € 
■p, Al > ... > A; > 0, we let rfc(A) = {j : Xj = k} to denote the number of parts 
of A that are equal to A: > 1. Then /i„ is defined by 

M.({Ae7':r,(A)=j})= ('^''^ + ^ " ^')^;'=J"(1 - ^;'=)''^ < « < 1. 

The key feature in the study of this kind of measures is the fact that a kind 
of a conditional probability measure on Vb{n) turns out to be independent of 



V for all n and coincides with the uniform probability measure P = Pn,b (for 
more details, see [19]). In [20] Vershik and Yakubovich studied the limiting 
distribution of the maximal particle energy X{X), or, which is the same, the 
largest part size X(A),A e V, with respect to the measure /z„, as w — > 1~. In 
particular, under the assumption that the weight sequence b satisfies (|1.5p . they 
proved that 

lim /i„({A eV :{1- v)X{X) -iy\log{l~v)\ 
-{u- l)log| log(l-t;) I - {i^-l)\ogi^-\ogC <t}) 
— e^'^ , — oo < t < cxo. (1-6) 

As it was mentioned before, the weight sequence (|1.5|) satisfies Meinardus con- 
ditions (Ml) — (M3). Vershik and Yakubovich [20] studied also a more realistic 
model of quantum ideal gas in a d-dimensional space, for which the weights 
satisfy J2i=i ^i — Cdk"^^^ + 0{k^''), as fc — > cxo {cd and kj, < d/2 are computable 
constants) . 

The main result of this paper is obtained in terms of the uniform probability 
measure P = P„^f, on the set Vb{n)- Before stating it, for the sake of brevity, 
we introduce the following notation: 

a„(p,A)= , n > 1, (1.7) 



where the constants p and A are defined by condition (Af 1). 

Theorem 1 If the weight sequence b, satisfies conditions {Ml) — (Af3), then, 
for all real t, the limiting distribution of the largest part size Xn is given by 



hm P{anXn + p log a„ - (p - 1) log | log a„ | - (p - 1) log p-logA<t) 

(1.8) 



^ — e 



Remark. One can easily compare (jl.81) with (jl.6p setting in the latter one 
V = 1 — QmV = p and C = A and observe the coinciding normalizations. We 
also notice that the limiting results for linear and plane partitions (see [6,12]) 
follow from (|1.8|) with A — 1 and p = 1 and 2, respectively. 

The method of our proof combines Hayman's theorem for estimating coeffi- 
cients of admissible power series [10] (see also [7; Sect. VIII. 5], a generalization 
of Perron's formula, expressing partial sums of a Dirichlet series by a complex 
integral of the inverse Mellin transform applied to the Dirichlet series itself (see 
Thm 3.1 from the Supplement of [16]) and some Mellin transform computations. 

We organize our paper as follows. Section 2 includes some auxiliary facts 
that we need further. Some proofs are omitted since they are given in [9,11]. 
In Section 3 we present the proof of Theorem 1. The Appendix contains some 
technical details related to the application of the generalized Perron's formula 
[16]. 



2 Preliminary Results 

We start with a lemma establishing an asymptotic estimate for infinite product 
representation (jl.f I) of the generating function fb{x). It has been proved by 
Meinardus [11] (see also [2; Lemma 6.1]). 

Lemma 1 Suppose that sequence b is such that the associated Dirichlet series 
Ijl.S]) satisfies conditions (Mi) and {M2). If t = a -\- iO , then 

fbie-n = exp {AT{p)C{p + 1)t-p - D{0) logr + 0(a^«)) 

as a — >■ 0+ uniformly for \ \< t: and \ argr |< 7r/4. ■ 

Our first goal is to show that Meinardus conditions (Mi) — (M^) imply that 
the generating function fb{x) possesses Hayman's admissibility properties [10] 
(see also [7; Sect. VIII.5]) in the unit disc. Hence, for < r < 1, we introduce 
the functions: 

oo 

Fbir)=\ogfbir) = -J2bk\ogil^r''), (2.1) 



fc=i 



Jiir) 



Ab{r)^rFl{r)^rf^ (2.2) 

fb{r) 

Bb{r) = r'Fl'ir) + rFl[r) =rlM+ r^f^ - r^ i^-j^S ■ (2-3) 

Ib[r) fb{r) \fb[r)J 

Furthermore, setting in (|2.ip - (|2.3p r = e~", we shall obtain their asymptotic 
expansions as a — ^ 0+. For the sake of convenience, we also set 

h = hip,A)^AT{p+l)C{p + l). (2.4) 

The proof of the next lemma is contained in [9; Lemma 2]. 

Lemma 2 Meinardus conditions (Mi) and {M2) imply the following asymptotic 
expansions: 

Abie-") = ha-P-^ + D{0)a-^ + Oia^°-^), (2.5) 

Bb{e-") = -^(-A(e-")) = h{p + l)a-P-^ - D{Q)a-^ + 0(a^°-2)^ ^3.6) 
da 

Fl"^Oia-P-^), (2.7) 

as a ^ 0"'", where Fb,Ab and Bb are defined by Ii2.1\) - ^27B) . respectively. More- 
over, from \2. 5\) it follows that the equation 

Ab{e-°')^n, n>l, (2.8) 

has a unique solution a — an, such that a„ —^ as n —>■ cxi. An asymptotic 
expansion of this solution, as n -> 00, is given by 

an ^ a„ + , •^^°^. +0{n-^-^), (2.9) 

(p+ l)n 

where /3 = min I -^?j-, ^^ ) and a„ are the normalizing constants given by |J.7|). 



We notice that ^^, (HSl) and (US]) imply that 

A(e~"")-^oo, 66(e^""^oo, n^oo, (2.10) 

that is, Hayman's "capture" condition [7; p. 565] is satisfied with r — rn — e~"" . 
Our next step is to establish Hayman's "locality" condition, which implies the 
asymptotic behavior of fb{x) in a suitable neighborhood of a; = 1. 

Lemma 3 Suppose that the weight sequence b satisfies Meinardus conditions 
(Ml) and {M2) and a„ is the solution of \2. i^|) given by 112. 9\) . Let 

5n = al+P/'^/uin), n > 1, (2.11) 

where uj{n) — > 00 as ri —> 00 arbitrarily slowly. Then 



(2.12) 



(2.13) 



^-Wn h{e °"+'' ) ^ ^_e^s.(e--.)/2(i + 0(l/c,3(„)) 

/b(e "") 

uniformly for | |< (5„. 

Proof. Applying Lemma 1, we observe that 

-.,„ A(e-""+-^) 
/fc(e-"") 

= exp ( -((a„ - iff)-" - a-P) - D{Q) log ( 1 - — ) - iOn + 0{a^° 

where h is given by p.4p . Expanding {an — i9)^P and log (1 — id jotn) by Taylor's 
formula and using ()2.5p . p.6p and p.Sp . we obtain 

-((a„ - i6')-'' - off) - ^(O) log f 1 - — V «^" 
p V an/ 

= ^e{ha-r^ + i?(0)«-i - n) - ^/.(p + X)^;/-^ - ^^^ + 0(| 1^ a.T^-") 
= i0(A(e-"") - n + O(a^o-i)) - y (6,(6—) + 0(a^''-2)) + 0(| 1^ a-^^") 

= -^6,(e-"") + 0(5„a^"-i) + 0(<5^a^"-2) + 0{bla-J-P). 

Substituting this into p.l3p and taking into account p. lip , we obtain (|2.12p . ■ 
To complete our analysis, we also need to study the behavior of /&(e~"""''**) 
outside the range —8n < 9 < 6n. The next lemma shows that Hayman's last 
("decay") condition [7; p. 565] is also valid. 

Lemma 4 Suppose that fb{x) satisfies Meinardus conditions (Afi) — (Afa). Then, 
for sufficiently large n, 



|/,(e-""+^^)|</fc(e-"")e- 



-Caa" 



uniformly for Sn ^ & < tt, where C2 , e > are the constants defined in condition 
(Ms). 
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Proof. First, we notice that 

exp(5R(log/6(e-""+"')) -log/,(e-"")). (2.14) 



/.(e-""+^^)| _,^„^„,,„_„„+,:e 



A(e-"") 

Then, setting 9 = 2ttu, we almost repeat the argument from [9;p. 324 
3?(log/,(e-""+''^))~log/,(e-"") 

/ °° /I „-a„+27rm\ \ 



■ 2^ fefc log 



Iv^, , / 1 - 2e-*"" COS (27ru/c) + e^ 



1 v^ , , / 4e-""''~ sin^ (Trufc) 



fe=i ^ ^ ' 

-. CO 

< -- ^ 6fe log (1 + 4e-""'= sin^ (Trufc)) 

fc=i 

< -^ X] ^fee-""'^ sin^ inuk), (2.15) 

fc=i 

where the last inequality follows from the fact that log {I + v) > ( -^^ ) v,0 < 
V < 4. Substituting this estimate into (I2.14p and applying condition (Mg), we 
obtain the required inequality. ■ 

We now recall (|2.6p from Lemma 2. It implies that 

Combining this asymptotic equivalence with the result of Lemma 4, we obtain 
Hayman's "decay" condition [7; p. 565], namely, 

I A(e"""+*^) h o(/,(e"")/Sy'(e-"")), ^ ^ oo, (2.16) 

uniformly for (5„ <| |< tt. 

Eqs. (|2.10p . (|2.12l) and (j2.16p show that the function fb{x) is admissible 
in the sense of Hayman. Therefore, we can apply Thm. VIII.4 of [7] for its 
coefficients. We state this result in the next lemma. 

Lemma 5 Suppose that the weight sequence b satisfies Meinardus conditions 
(Ml) — (M3). Then, the asymptotic for the total number of weighted partitions 
is given by 

^ e- """/,(e-"") 

v/27rS,(e-"") 



as n ^ 00, where a„ is the unique solution of 112. 8\) whose asymptotic expansion 
is given by 112. 9]) and Bb{e~"") is defined by 117 



Remark. The asymptotic equivalence (I2.17P is in fact Meinardus asymptotic 
formula [11] for the number of weighted partitions of n. It was also established 
by Granovsky et al. [9] under a condition weaker than (M3). Here we give the 
formula in a slightly different form, which is more convenient for our further 
asymptotic analysis. One can easily show the coincidence of (I2.17P with the 
Meinardus original formula, applying the result of Lemma 1 to /^(e^"") and 
replacing a„ and ;Bb(e~"") by (j2.9p and (|2.6p . respectively. 

Further, we also need a bound on the rate of growth of the weights b^, as 
fc — > 00. Granovsky et al. [9; p. 310] obtained it using Tauberian theorem 
technique. 

Lemma 6 If the sequence of weights b satisfies Meinardus conditions (Mi) and 
{M2), then 

bk = o{kP) 

as fc —T' 00. ■ 

Now we recall formula (|1.2p for the truncated products fm.b{x),m > 1. 
Similarly to (|2.ip . we set 

m 

F,r,,b{x) = log/,„,b(a;) = -^6felog(l - x*'). (2.18) 

fc=i 

(Here we consider the main branch of the logarithmic function, assuming that 
logy < for < y < 1). Further on, when computing the derivatives of (j2.ip 
and (|2.18p . we shall write 

F(^")(e-"") = fI^\x) !_,-.„, Fi^;i(e--) = F^ (x) U.e-=.,J = 1,2,3. 

Our next lemma establishes estimates on the tails Fjf'{e^°'") — -?„ b(e~"") for 
some specific values of m. 

Lemma 7 Suppose that the weight sequence b satisfies Meinardus conditions 
(Ml) and (A/2) and a„,n > 1, is defined by eq. 112. 9\) . Moreover, let 

m ^ pa~^ log a~^ , n — >• 00. (2.19) 

Then 

F(^-)(e-.) _ /^«,(e-"") = 0{a-^ log^+^' a.-^), j = 1,2,3. 

Proof. We shall consider the case j = 1. The other two cases are studied in 
a similar way. 

First, we choose an integer mi that satisfies the asymptotic equivalence 

mi^(p+l)a,;Moga-i (2.20) 



n'(e-"")-^:,.(e-"")= E 1 ,-^".^ -^1+^2, (2.21) 



and decompose the difference of the first derivatives in the following way 

k—m-\-l 

where 



p — r^uLji >• -< I p — ko 

/c— m+1 k—rui +1 



From Lemma 6, (|2?T9| and (|2?20l) it follows that 

1 '"1 -(m+l)Q„ 

= o(mf ^e-™"") = o(a-(''+i)(log''+i a-^)a(:) 
^oCa-Mog^+^a-i). (2.22) 

S'2 can be estimated in a similar way, using the asymptotic expansion of the 
incomplete gamma function r(a, z) as z —> 00 (see [1; Sect. 6.5]). We have 

= O (a-P-^ I uP+^e-'^dv\ = 0(a-''-2(mia„)''+ie-'"i"") 

= 0(a-ir<ie-(''+i'i°8""') = 0(a-Mog''+i a-^). (2.23) 

The required estimate now follows from (|2.21I) - (|2.23I) .B 

Our last forthcoming lemma supplies us with integral representations for 
F(e~") and i^m(e~"),a > 0, using Dirichlet series p.Sp and its partial sums 

771 

Dra{s) ^^hkk-\ s^a + iy, m > 1, (2.24) 

respectively. The proof is based on Mellin transforms and can be found in [11], 
[9; Lemma 2(ii)] and [2; Sect. 6.2]. 



Lemma 8 For any a, A > 0, we have 

/•p+A+ioo 

/p + A — 700 

and 



1 i-p+A+ioo 

FmAe~n = TT- / a-'T{s)Cis + l)D„^{s)ds (2.25) 



1 pp-i-i^+ioo 

Fb{e-") - :— / a-'r{s)as + l)D{s)ds, (2.26) 

27ri Jp^A-ioo 

where Dm{s) and D{s) are defined by \2.24^ and hl.S^) . respectively. M 
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3 Proof of the Main Result 

We apply first Cauchy coefficient formula to (jl.2p using the circle x = e""""*"**, tt < 
< TT, as a contour of integration (a„ is determined by (|2.9p V We obtain 

Pb{n)P{X„ < m) - —— / /„,,b(e-""+^«)e-*«"d0. 

Then, we break up the range of integration as follows: 

Pb{n)P{Xn <m) = Ji(m, n) + J2(m,n), (3.1) 

where 






2^ J-s^ /m,;,(e-"") 



"-de, (3.2) 



J.(m, n) = ^ / Y r-^ A "'"^^' (3.3) 

where (5„ is the sequence defined by (|2.1ip . 

We start with an estimate for Ji{m, n), expanding the integrand of p.2p by 
Taylor's formula: 

, ,^-^ . ' = exp{(e^^ - l)e-""F:.,(e--) 

Jm,b\S "■ ) 

+ i(e^« - l)2e-2""i^l,(e-"") + 0{\ 6 p F-,(e-""))}. 
Hence, we can rewrite (13.21) as follows: 



where 



Ji(m,n) = ^^e^-''(^"°")/„, (3.4) 

V 27r 



In = ^r exp{(e*« - l)e'""^:,,(e-"") 
V27r j-5„ 

+ 1 (e'« - l)2e-2«"F^ „(e-"") + 0{\ 9 \^ FZ,{e~'^-)) - ^en}de. 

Lemma 7 shows that, for those m satisfying (|2.19p . we can replace the derivatives 

-^m fc('^ "") ^y -Pfc (e~""), j = 1, 2, 3, at the expense of a negligible error term. 
In fact, combining Lemma 7 with (j2.1ip . we have 

(e^' - ^YF^lie-"") = ie'' - l)^F(^-)(e-"") + 0{5ia~^ log''+^- a'') 
= (e^" - l)^F,(^)(e-"") + 0(a«/3 log''+^' «,;V..^"(«)), .7 = 1, 2, 
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and 

0{\ e r <.(e-"")) = OK(log''+^a-i)/c.3(r.)), 

for any function w(n) — >■ cx) as n — >■ cxd arbitrarily slowly. Since all error terms 
above tend to 0, we obtain 

/„ = i±|^ /'" expKe^'' - l)e-""^,'(e-"") 



l + o(l) 



V27r J_5, 



/ exp{Ff,(e-""+''') - -Fb(e-"") - ien}de 



-i&n 



27r y_5„ /6(e-"") 



d6l. 



J-5„ 


1 




1 .5„Ve.(e-=") 


/•OO 


1 

T7 — 5. nn 


v/27r66(e- 


-"")y-oo 



where the last equality follows from a similar Taylor's expansion for F\,[e ""+'®). 
Now, from Lemma 3 it follows that 

In 



The last asymptotic equivalence follows from p.6p of Lemma 2, which implies 
that 

-p/6 

if w(n) — >■ oo slower than «„'' . 

Substituting the asymptotic equivalence for /„ into p.4p . we conclude that 

p-ncin 

whenever m satisfies ()2.19p as n — >■ oo. 

For the estimate of J2(m,n), we recall p.3p and the proof of Lemma 4. 
Using an argument, similar to that presented in p.lSp . for any real u, we obtain 

1 '^ ."'. 

3?(i^m,6(e-""+2™)) - F™,b(e-"") < -^ ^ fofce-""'^ sin^ {^uk) 

fe=i 

= -^ 5] 6fce~""'= sin^ (TTuft) - ^ fofee^""'^ sin^ (^ufc) . (3.6) 

\fe=l k=m+l J 
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Combining Lemma 6 with (|2.19p and repeating the argument given in the proof 
of Lemma 7, we also get the following estimate: 

oo oo 

Y, bkc-^'-'' sin^ (nku) < J2 &fce-""'= = 0(log''a-i). 

k—rn~\-l /c— m+l 

Replacing the second term of the right-hand side of p.6p by the last O-estimate 
and applying condition condition (Mg) to its first term, for 5„/27r <| u \< 1/2, 
we obtain 



Now, we are ready to compare the growth of (|3.3p with that of (j33|). So, if m 
satisfies (|2.19p as n — > cx), we have 

I J2im,n) |< exp(-na„ + F„,fc(e""")) 

X / I /™,6(e-""+""")//™,6(e-"") I du 

= exp(-na„ + i^„,b(e-"")) / (3?(i^™,fc(e-""+2-") - F,„,fc(e-""))dw 



= 0(e-'^^""V27r66(e~"")Ji(m,ri)) = o(Ji(m,n)). (3.7) 

It is now clear that p.ip . p. 51) and p. 71) imply that 

Pb{n)P{X„ < m) pP^Ae °")^ n ^ oo. 

^27r66(e-"") 

Subsequent application of the asymptotic equivalence (|2.17p from Lemma 5 
implies that 

PiXn <m)^ exp{i^„,b(e-"") - Fb{e'"-)}, (3.8) 

where q;„ and m satisfy (|2.9p and (|2.19p as n — ?> oo, respectively. 

Further on we shall study the asymptotic behavior of the exponent in p.Sp . 
Our analysis will be based on a generalization of Perron's formula that expresses 
partial sums of a Dirichlet series as complex integrals of the inverse Mellin type 
transforms applied to the Dirichlet series itself. We shall use it in the form 
given in the Supplement of [16; Sect. 3]. So, first we represent -F'm,b(e~"") using 
representation (|2.25p of Lemma 8 and then, we apply Perron's formula to the 
partial sum Dm{s) of the Dirichlet series D{s) that is included in the integrand 
of (|2.25p (recall also (|2.24p and (|1.3p '). In this way we arrive at the following 
complex integral representation: for any A > 1, we have 

F,n,b{e-^") (3.9) 

= — / a-^T{s)as + 1) ^ ^ ^ + D{s) + o(l) ds. 
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Furthermore, p.8p and p.9p imply that 

P{Xn < m) (3.10) 

~exp -— ^ / ((™+l)a„)-r(5 + p)C(s + p+l)- . 

[^ 27^^ 7a-joo s J 

The proofs of p.9p and p. 101) contain some technical details that will be given 
in the Appendix. 

We continue with the computation of the complex integral in the exponent 
of (|3.10l) . The exact value of m = m(n) will be specified later. First, setting in 
the integral of ((3?T0)) 

u = Un = (to + l)a„, (3.11) 

we consider it as inverse Mellin transform of the function T{s + p)C,{s + p+l)/ s. 
For the sake of convenience, we also set 

1 /-A+Joo , 

H{u) = — u'^ns + p)C{s + p+l)-. (3.12) 

27^^ Ja-joo s 

Then, clearly gi{s) — 1/s is the Mellin transform of 

„ , . { I if u < 1, 

^ ' \ if u < 1 

(see [5; formula 6.2.18]), while 32(5) — T{s)Q{s) is the Mellin transform of 



i?2(«) = £^— = -log(l-e-") 



.=1 ^ 



(see [7; p. 764]). Next, for A > 1, we apply formula (6.1.14) from [5] with a = 
and /? = p — 1. We obtain 

/■OO /"OO 

H{u) = u'^ / ^Hiiu/OH^iOd^ = - / y"-' log (1 - e~y)dy 

Jo Ju 

yP'^e~ydy + R{u) = r{p, u) + R{u), (3.13) 



where T{p,u) denotes the incomplete gamma function, while R(u) is an error 
term given by 



Riu) = j^ yP-' ^ __ dy. 
It is easily estimated as follows: 

/•OO 

R{u) = 0{ yP-^e-^ydy) = 0(e-"r(p, u)), u^ 00. 

Ju 
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Combining this estimate with p.l0p - p.l3p and applying the asymptotic expan- 
sion of the incomplete gamma function [1; Sect. 6.5], for u = w„ ~ man -^ oo, 
we obtain 

P(X„ < m) = exp{~Aa-P{uP-^e-"" + 0«-ie-2-"))} 

= exp{-^a-im''-ie-'"""(l + 0(e""""))}. (3.14) 

It is now clear that P{Xn < rn) converges to the distribution function e^"^ , — oo ■ 
i < cx), if TO = 771 (n) satisfies 

-man + (p - 1) log TO + log (Aan^) = -t + o(l) 

as n — >■ oo. From this we deduce 

TO = a~"^loga^^ + (p- l)a~"^logTO+ (logA + t)a~^ + o{a~^), (3.15) 

which in turn implies that 

logTO = log(a;;Moga;;i) 

/ log A + 1 , , log m 

+ log 1+ , _i +(p-l)-^- 

V log a„ log an 

= log a^ ^ -flog log a;^^ 

+ log('l+^"g-'V + (p-l)^"g""+^"g^"g^"+^(^ 

V log a„ log a„ 

= loga;:! +logloga-i +log fl + (p- 1) + O ("iSil^i^ 

1 -1,11 -1,1 , ^ /logloga,T^~ 
= loga„ -|-logloga„ -l-logp-fO — ^j- 

V log a„ 
Hence, p.l5p becomes 

TO = pa;;;^ log a;;^ -i- (p - i)^;;;^ log log a''^ 

+a-\p -I) log p+ {log A + t)a-'+o(^^f^^^ 

\ log an 

Replacing now this value of to into p.l4p and using the continuity of the dis- 
tribution function e^'^ , — oo < t < oo, we obtain 

P{Xn < to) 

= P{anXn - plogan^ - (p - l)logloga,7^ - (p- l)logp-log74 < t) + o{l) 
-^ e^"'', n^oo. (3.16) 

To complete the proof of the limit theorem (11.81) . it remains to justify the ap- 
propriate normalization replacing the sequence a„ , n > 1 , in p.l6p by a„ = 
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an{p,A),n > 1, given by (|1.7p . To show this we recall (|2.9I) and notice that 
after taking logarithms from its both sides, we easily obtain 

loga„ =loga„ + 0(n"^). (3.17) 

Furthermore, (|2.9|) implies that 

O-nXn — anXn — i^ri , (3.18) 

where ¥„ — , !^J^ Xn- From (|3.16p it is not difficult to conclude that, for 
every 77 > 0, lim„_>oc P{Yn > i]) ~ 0, that is, the sequence of the random 
variables Yn,n > 1, tends to in probability, as n — ^ 00. In other words, 
p.l8|) represents a„X„ as a sum of two random variables: the first one, anXn-, 
converges in distribution, as n — )> cx), to a random variable with a distribution 
function e~'^ ,— oo<t<oo, while the second one, Y„, tends to in probability. 
Using this fact, (|3.17p and the continuity of the distribution function e"*^ , it 
is an easy probabilistic exercise to show that the limit theorem given by (|1.8p 
holds. 



4 Appendix 

Proof of i3.9]) . First, we recall eq. (|2.25p of Lemma 8. Our goal is to represent 
the TTith partial sum Dm{s), defined by ()2.24|) . using the inversion formula given 
by Thm. 3.1 in the Supplement of [16] (see also formula (3.4) there). Instead 
of D{s), we shall consider now the Dirichelet series 

00 

k=l 

which converges absolutely for ct > 1. Lemma 6 implies that the coefficients of 

D{s + p — 1) satisfy 

bkk-P+^ = o{kP)k-P+^ = o{k) < dk 

for any constant c > and all k (in other words, the function $(a;), introduced 
in [16] satisfies $(x) = a; in this case). Furthermore, from Meinardus condition 
(Ml) it follows that 

^ A 

J2bkk-P^'k- 



cr- 1 
fc=l 



where (f>{s) denotes a function which is analytic for a > —Co- Hence 
00 
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So, the conditions of Thni. 3.1 (the Supplement of [16]) are satisfied and by 
its second part we conclude that, for large enough T > and any pair of fixed 
constants A > 1 and d > 0, we have 

D„^{w + p-l) + ^b„,+i{m + l)-P+\m + l)""' 
= / D(w^z^p-\)- -dz 

+0 I — ] + O I — -^ — \ ,w = 1 + A + iv, -oo < w < oo. (4.1) 

Lemma 6 implies that the second term in the left-hand side of (|4.ip is o{m^^P^^) 
as m — ^ cx). To compute the integral in the right-hand side of (14. ip . we use a 
contour integral around the rectangle d — iT, d + iT, — Co — p — A + iT, —Co—p — 
A — iT. Using Meinardus condition (M2), we estimate the integral over the end 
segment (-Co - p- A + iT, -Co - p- A-iT)hy 0{T^'m-^°-P-^). Hence, 
it tends to as ttt., T — > cxj, provided 

r = o(m(^°+''+^)/cri)- (4.2) 

An estimate for the integrals on the segments (— Cq — p — A + iT, d + iT) and 
(—Co — p — A — iT, d — iT) is given by 

, m'^ \ /to' 



^TlogmJ \ T 

Taking into account the first 0-estimate in the right-hand side of (|4.ip . we 
conclude that T and d should also satisfy 

m'^^o{T), m,T^ 00, (4.3) 

Therefore, further we shall assume that to, and T satisfy (|4.2I) and (14. 3p with 
certain constants A > 1 and d > 0. Thus, all integrals on the end segments 
except the integral on (d — iT, d + iT) are negligible for such to and T as ?7i, T — >■ 
00. The same obviously holds for both O-estimates in the right-hand side of 
(j4.ll) . So, we can compute D„i{w + p — 1) summing up the residues of the 
integrand in (|4.ip . We obtain 

Aim -\- 1)"'^^™ 

D„i{w+p-l) = — \-D{w+p-l)+o{l),w = l+A+iv,-oo < u < 00. 

1 — w 

Setting w — s — p + 1 and substituting this expression into p.25p . we arrive at 
dSH). 

Proof of ( TOO)) . dSH) and ([2:261) imply that 

F™(e-"") = — / a-^r(s)C(s -f 1)^^ ds 



27rz Jp+A-ioo 



+F(e-"") + o(l) 

((to + l)an)-'T{s + p)C{s + p+l)—+ F(e-"") + o(l). 



Aa-P f^+'°° __, _, ^ds 



27ri Ja-xoo s 
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Replacing this into p.Sp . we obtain p.lOp . 
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